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The magnetoresistance of a granular superconductor in a strong magnetic field is considered. It 

O ' is assumed that this field destroys the superconducting gap in each grain, such that all interesting 

C~^ [ effects considered in the paper are due to superconducting fluctuations. The conductance of the 

system is assumed to be large, which allows us to neglect all localization effects as well as the 

Coulomb interaction. It is shown that at low temperatures the superconducting fluctuations reduce 

the one-particle density of states but do not contribute to transport. As a result, the resistivity 

of the normal state exceeds the classical resistivity approaching the latter only in the limit of 

C^ ' extremely strong magnetic fields, and this leads to a negative magnetoresistance. We present detailed 

i-i^ , calculations of physical quantities relevant for describing the effect and make a comparison with 

iy>j ' existing experiments. 



^ 



CD 

B 



o 



a 



^ 

a 



PACS numbers: 73.23.-b, 74.80.Bj, 74.40. +k, 72.f5.Rn 



I. INTRODUCTION 



Ph ' In a recent experimentEl, transport properties a system of Al superconducting grains in a strong magnetic field were 

O ' studied. The samples were quite homogeneous with a typical diameter of the grains 120 ± 20A and the grains formed 
a 3-dimensional array. As usual, sufficiently strong magnetic fields destroyed the superconductivity in the samples 
and a finite resistivity could be seen above a critical magnetic field. The applied magnetic fields reached 17T, which 
^~^ • was more than sufficient to destroy also the superconducting gap in each grain. 

^ ■ The dependence of the resistivity on the magnetic field observed in Ref. EJ was not simple. Although at extremely 

^.. strong fields the resistivity was almost independent of the field, it increased when decreasing the magnetic field. Only 

i«— » ' at sufficiently weak magnetic fields the resistivity started to decrease and finall}j-|tjl|ie samples displayed superconducting 
f— ^ , properties. A similar behavior had been reported in a number of publicationaS'Ll. p. 

y—i A negative magnetoresistance due to weak localization effects is not unusual in disordered metalaj. However, the 

0^ , magnetoresistance of the granulated materials studied in Rcf.El is quite noticeable in magnetic ffelds exceeding lOT 
On ■ and is considerably larger than values estimated for the weak localization. 

The weak localization effects become very important if the system is near the Anderson metal-insulator transition 

and one can expect there a complicated dependence on the magnetic field. One can, for example, arguea that decreasing 

the magnetic field drives the system to the metal-insulator transition, which gives a large negative magnetoresistance. 

If the superconductimK. transition occurs earlier one gets a maximum of the resistivity near the transition characteristic 

for the experimentsEI tl. 

Q Unfortunately, investigation of this possibility is not simple. The metal-insulator transition occurs at values of the 

O ■ macroscopic conductance J of the order of unity. At such values calculations arc very difficult. The problem becomes 

^ I even more complicated due to the Coulomb interaction. It is well knowriJ, that, at small values of J, the system 

• ^H . must be an insulator even if the superconducting gap is finite in a single grain. A microscopic consideration of all 

^\ ' these effects and a confirmation of the existence of the negative magnetoresistance near the superconducting point 

'^ for J ~ 1 is at present hardly possible and we do not try to treat the problem here. 

■ - - ' Instead, we consider below the region of large conductances J ^ 1, where the system without interactions would 

be a good metal. This region corresponds to large tunneling amplitudes between the grains. All effects of the weak 

localization and the charging effects have to be small for J :s> 1, which would imply that the resistivity could not 

considerably depend on the magnetic field. 

Nevertheless, we find that the magnetoresistance of a good granulated metal (J ^ 1) in a strong magnetic field 
and at low temperature must be negative. In our model, the superconducting gap in each granule is assumed to be 
suppressed by the strong magnetic field. All the interesting behavior considered below originates from the supercon- 
ducting fluctuations that lead to a suppression of the density of states (DOS) but do not help to carry an electric 
current. The main results have been already publishedQ and now we want to present details of calculations and 



clarify some additional questions. We consider, for example, influence of the Zeeman splitting on the resistivity, flnd 
the critical field Hc2 of the transition into the superconducting state, calculate the diamagnetic susceptibility and 
estimate weak localization corrections. 

Theory of superconducting fluctuations near the transition into the superconducting state has been developed long 
agoLTu (for a review see Ref. nffl). Above the transition temperature Tc, nmji-cquilibrium Cooper pairs are formed and 
a new channel of charge transfer opens (Aslamazov-Larkin contribution)Ll. Another fluctuation contribution comas 
from a coherent scattering of the electrons forming a Cooper pair on impurities (Maki-Thompson contribution)El. 
Both the fluctuation corrections increase the conductivity and, when a magnetic fleld is applied, lead to a positive 
magnetoresistance. Formation of the non-equilibrium Cooper pairs results also in a fluctuational gap in the one- 
electron spectrumEl but in conventional (non granular) superconductors the flrst two mechanisms are more important 
near T^ and the conductivity increases when approaching the transition. The total conductivity for a bulk sample 
above the transition temperature Tc can be written in the following form 

a = (TDrude + 5a, (1.1) 

where (TDrude — {e'^Tn)/m is the conductivity of a normal metal without electron-electron interaction, r is t he elastic 
mean free time, m and n are the effective mass and the density of electrons, respectively. In Eq. ( p, . l| ) , Sa is a 
correction to the conductivity due to the fluctuations of the virtual cooper pairs 

Sa = Sanos + SaAL + SaMT (1-2) 

where Saoos is the correction to the conductivity due to the reduction of the DOS and 6a al and 5a mt stand for the 
Aslamazov-Larkin (AL) and Maki-Thompson (MT) contributions the conductivity. Close to the critical temperature 
Tc the AL correction is more important than both the MT and DOS corrections types and its contribution can be 
written as followsQ 

^-^^ x(-I^X (1-3) 



where A is a small dimensi onle ss positive parameter A <C 1 and f3 = 1/2 for the three-dimensional case (3D), 1 for 2D 



and 3/2 for quasi-lD. Eq. (1.3) was derived using a perturbation theory and therefore is valid provided the inequality 
SaAL/o-Drude < 1 is fulfilled. 

Although typically the AL and MT corrections are_larger than the DOS contribution, a small decrease of the 
transverse conductivity is possible in layered materialsliil in a temperature interval not very close to the transition. 
It is relevant to emphasize that all previous study of the fluctuations has been done near the critical temperature 
Tc in a zero or a weak magnetic fleld. In contrast, we are mainly concentrated on study of the transport at very 
low temperature in a strong magnetic fleld. To the best of our knowledge, fluctuations in this region have not been 
considered so far. 

A strong magnetic fleld destroys the superconducting gap in each granule. However, even at magnetic flelds H 
exceeding the critical fleld He virtual Cooper pairs can still be formed. It turns out, and it will be shown below, the 
influence of these pairs on the macroscopical transport is drastically different from that near Tc- The existence of the 
virtual pairs leads to a reduction of the DOS but, in the limit T — > 0, these pairs cannot travel from one granule to 
another. As a result, the conductivity a can be at i? > He considerably lower than conductivity ctq of the normal 
metal without an electron-electron interaction. It approaches the value ctq only in the limit H 3> H^ when all the 
superconducting fluctuations are completely suppressed by the magnetic field. 

The superconducting pairing inside the grains is destroyed by both the orbital mechanism and the Zeeman spliting. 
The critical magnetic field H°^ destroying the superconductivity in a single grain in this case can be estimated as 
H°^R£_ « (/)o, where (/)o = hc/e is a flux quantum, i? is a radius of single grain and f = \/^ is the superconducting 
coherence length. The Zeeman critical magnetic field H^ can be written as gfiBH^ = Aq, where Aq is a BCS gap for 
the single grain at magnetic fleld H = and g is a Lande factor. We notice here that H^ is independent of the size 
of the grain whereas for H°^ the size of the grain is important. The ratio of this two flelds can be written in the form 

H^'^/H^ « Rc/R, (1.4) 

where Re = C(PoO~^- We can see from Eq. ( |1.4| ) that for R > Re the orbital critical magnetic fleld is smaller than the 
Zeeman critical magnetic fleld H"^ < He and the suppression of sufjerconductivity is due to the orbital mechanism. 
This condition is. well satisfled in grains with R ~ lOOA studied inld. This limit is opposite to the one considered 
recently in RefJi3 where the Zeeman splitting was assumed to be the main mechanism of destruction of the Cooper 
pairs. However, the latter mechanism of the destruction of the Cooper pairs can be easily included into the scheme 
of our calculations. 



The remainder of the paper is organized as follows. In Sec. ^we formulate the model and discuss the fluctuational 
contributions to the total conductivity of the granulated superconductors. Sec. Ill contains the derivation of the 



correction to the conductivity of granulated superconductors due to single electron tunnelling (DOS contribution) 



at very low temperatures T <^ Tc and strong magnetic fields H — He ^ He- In Sec. IV and Sec. ^ corrections 
to the conductivity due to tunneling of virtual Cooper pairs are derived (Aslamazov-Larkin and Maki-Thompson 
corrections) . In Sec. VI we discuss the influence of magnetic field on the phase of order parameter and calculate the 
criti cal fi eld of the transition into the superconducting phase. The importance of the Zeeman splitting is discussed in 
Sec. |VIl|. T he contribution of fluctuating Cooper pairs to the diamagnetic susceptibility of granular system is derived 
in Sec. VIII. Sec. IX includes the results for conductivity at T — Tc <C Tc and H ^ He- A discussion of a recent 
experiment in Al grains and a comparison with the theory is presented in Sec. ^ Our results are summarized in the 
Conclusion. 
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FIG. 1. System of metallic grains 



II. CHOICE OF THE MODEL 



We consider a 3D array of superconducting grains coupled to each other, Fig. n^. The grains are not perfect and 
there can be impurities inside the grains as well as on the surface. We assume that electrons can hop from grain to 
grain and can interact with phonons. 

The Hamiltonian H of the system can be written as 



H = Hu 



Ht, 



(2.1) 



where Hq is a conventional Hamiltonian for a single grain with an electron-phonon interaction in the presence of a 
strong magnetic field 



Hq = 2_^^i,kai,k'^i,k — |A| 2_^ 0,i_kO,i_-k'^i,~k'ai,k' + Hi, 



(2.2) 



i.k 






where i stands for the numbers of the grains, k = (k, "f), —k = (— k, |); A is an interaction constant, and Himp 
describes elastic interaction of the electrons with impurities. The interaction in Eq. (2.2) contains diagonal matrix 
elements only. This form of the interaction can be used provided the superconducting gap Aq is not very large 



Ao<Sc 



(2.3) 



where Ec is the Thouless energy of the single granule. Eq. (2.3) is equivalent to the condition R ^ ^o, where R is 
the radius of the grain and ^o is the superconducting coherence length. In this limit, superconducting fiuctuations in 
a single grain are zero-dimensional. _ 

The term Ht in Eq.(2.1) describes tunneling from grain to grain and has the form (see e.g. Ref.Ej) 



Ht= 2_^ typ5aJpaj,exp(i-Ady) 



h.c. 



(2.4) 



where A is the external vector potential; d^j are the vectors connecting centers of two neighboring grains i and j 
(|dij| — 2R); a- (uip) are the creation (annihilation) operators for an electron the grain i and state p. 



It is assumed that the system is macroscopically a good metal and this corresponds to a sufficiently large tunneling 
energy t 



t>(5 



(2.5) 



where S — {voV) is the mean level spacing in a single granule, vq = mpo/2Tr'^ is the DOS per one spin in the 
absence of interactions, V is the volume of the granule, and Aq is th e Bardeen-Cooper-Schrieffer (BCS) gap at T = 
in the absence of a magnetic field. Provided the inequality (2.5) is fulfilled localization effects can be neglected 
Ref.tj. Moreover, charging effects are also not important in this limit because at such tunneling energies the Coulomb 
interaction is well screened. 

The tunneling current operator is 



j = -led ^ UjpqalpCLjq exp(i-Adij) - h 



(2.6) 



Using standard formulae of the linear response theory we can write the current j {t) in the form 



JW = * / m)r3{t')\)Ht')dt'~ie'd' / {[HT{tlHT{t')])dt'A{tl 



(2.7) 



where the angle brackets stand for averaging over both quantum states and impurities in the grains. All operators in 
the right hand side of Eq. (2.7) are independent of the vector potential. In principle, the grains can be clean and the 
electrons can scatter mainly on the surface of the grains. However, provided the shape of the grains corresponds to a 
classically chaotic motion of the electrons, the clean limit should be described in the QD case by the same formulae. 

We carry out the calculation of the conductivity making expansion both in fluctuation modes and in the tunneling 
term Ht- This implies that the tunneling energy t is not very large. Proper conditions will be written later but now 
we mention only that the tunneling energy t will be everywhere much smaller than the energy Ec- 

As in conventional bulk superconductors we can write corrections to the classical conductivity as a sum of corrections 
to the DOS and of Aslamazov-Larkin (AL) ual and Maki-Thompson (MT) umt corrections. Diagrams describing 
these contributions are represented in Fig. 2. 




FIG. 2. Diagram a) describes correction to DOS, diagrams b) and c) describe corrections to conductivity due to super- 
conducting fluctuations. The wavy lines denote the propagator of the fluctuations, the dashed lines stand for the impurity 
scattering. 



The total conductivity <t can be written as 

cr = (tdos + ^AL + cfmt 
where odos is given by equation 

+ 00 



CTDOS = ao (4r)-^ j He)/vof cosh-2(^)de, 



(2.8) 



(2.9) 



In Eq. (2.9), ctq = 27re^i? ^ [t/S) is the classical conductivity of the granular metal. It can be rewritten in terms of 
the dimensionless conductance of the system J as 



8Je2 
.0 ^ ^ (2.10) 

where the conductance J equals 

J^(7rV4)(i/<5)2 (2.11) 

The inequaUty ( p.q ) is equivalent to the condition 

J>1 (2.12) 



The function i/ {e) in Eq. ( |2.9D is the density of states. Without the electron-electron interaction, this function is 
equal to the DOS of the ideal electron gas i^o, which gives ctdos = <^o- Taking into account the electron-electron 
attraction we can write the contribution crjjos to the classical conductivity as 

aoos =o-o+SaDos{t,T,H) (2-13) 

The correction daoos (^: T, H) depends on temperature T and magnetic field H and is represented in Fig. 2a. As 
concerns the long range part of the Coulomb interaction (charging effects), the condition J 3> 1 allows us to neglect 
it. 



Using Eq. (2^) the correction to the conductivity Saoos at low temperatures can be written in terms of the 



correction to the DOS at zero energy Sv (0) as 

Saoos /'yo = 2(3:^ {0)/i^o) (2.14) 

As we will see below, the main contribution to the conductivity due to the superconducting fluctuations comes from 
the change of the DOS Sv. Its calculation will be presented in detail in the next Section. 

III. SUPPRESSION OF THE CONDUCTIVITY DUE TO DOS FLUCTUATIONS 



In this Section we consider the correction to the conductivity of granulated superconductors due to suppression of 
DOS. The main correction Sv (s) to the DOS of the non-interacting electrons vq is described by the diagram in the 
Fig. 2a, while the terms (Jmt and (tal are given by Figs. 2b and 2c, respectively. The calculation of the diagrams 
can be performed for the-^atsubara frequencies e„ = 7rT(2n -I- 1) using temperature Green functions. At the end 
one should, as usual Ref.li3, make the analytical continuation ie„ -^ e . The magnetic field will be considered in the 
quasi-classical approximation I <^ Lc, where I is a mean free path and Lc is a cyclotron radius. In this approximation, 
the magnetic field results in the appearance of additional phases in Green functions 



G{ien,r,r')=G^"\ien,r-f')e^p \- Adf\, (3.1) 



where G''*''(ie„, r — r') is the Green function without magnetic field. In the zero order approximation in the super- 
conducting fluctuations, the disorder averaged Green function G'"' in the momentum representation has the form: 

G(°) (i£„, p) = (isn - e (p) + i {2t)-' sgnsn) ' (3.2) 



The diagrams in Fig. 2 contain the averaged one-particle Green functions, the impurity vertices proportional to the 
so called Cooperon C and the propagator of the superconducting fluctuations K. The functions C and K depend 
on the coordinates and time slower than the averaged one-particle Green functions because the characteristic scale 
for both the impurity vertices and the propagator of superconducting fluctuations K is the coherence length ^ which 
is much larger than I. As a result, the magnetic field affects only the vertex C and the propagator K, whereas the 
phases of the Green functions drawn in Fig. 2 outside these blocks cancel. So, reading the diagrams in Fig.l one should 
replace the solid lines by the functions G^'^' {iSn, p). More complicated diagrams containing crossings of impurity lines 
describe the weak localization effects are neglected here. 

The impurity vertex entering these diagrams is equal to {2tti>qt) C (iSmi^k — *£«), Fig. 3, where r is the mean 
free time due to the scattering on impurities or on the grain boundary, C is the Cooperon. It obeys the following 
equation 



(^Do i-iV- (2e/c) A)' + |2e„ - Qk]) C (r, r') = 2^^.0-5 ( 



(3.3) 



where Dq — Vqt/3 is the classical diffusion coefficient. The vector-potential A (r) should be chosen in the London 
gauge. If the shape of the grain is close to spherical, the vector-potential is expressed through the magnetic field H 

as A(r) = [Hxr]/2. 
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FIG. 3. Impurity vertex. 



All relevant energies in the problem are assumed to be much smaller than the energy of the first harmonics 
Ec ~ DoTT^/R^ playing the role of the Thouless energy of a sing le g rain and this allows to keep only the zero 
harmonics in the spectral expansion of the solution C (r, r') of Eq. (3.3). One can find the eigenvalue Eq (H) of this 
harmonics using the first order of the standard perturbation theory 



£o{H) 



(2e/cf Do < A2 >o 



(3.4) 



where < 
obtains 



>o stands for the averaging over the volume of the grain. For the grain of a nearly spherical form one 



So (H) = 



eHRV 



Dn = 



TT^o 



E, 



where (J)q — nc/e is the fiux quantum and is the magnetic fiux through the granule. 

Within the zero-harnronics approximation, the function C does not depend on coordinates and equals 

C{ien,i^k-i£n) = 2TTVo{\2£n-^k\ +£o{H)y^ 



(3.5) 



(3.6) 



To calculate the propagator of the superconducting fiuctuations K one should sum the sequence of the ladder 
diagrams represented in Fig. 4. The broken lines in this figure denote the electron-electron interaction. As it has 
been mentioned the characteristic energies of the propagator K are low and therefore, when calculating the function 
K, o ne should take into account the tunneling processes from grain to grains. The tunneling Hamiltonian iJy, Eq. 
(2.4), is represented in Fig. 4 by crossed circles. 



K 
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FIG. 4. Propagator of superconducting fluctuations with tunneling. 



Of course, one can sum the ladder diagrams in Fig. 4 directly. However, sometimes it is more convenient to decouple 
the electron-electron interaction in Eq. (pj) by a Gaussian integration over an auxiliary field A (Hubbard-Stratonovich 
transformation). Then, one can perform averaging over the electron quantum states, thus reducing the calculation to 
computation of a functional integral over the field A. In principle, one obtains within such a scheme a complicated 
free energy functional and the integral cannot be calculated exactly. The situation simplifies if the fiuctuations are 
not very strong. Then, one can expand the free energy functional in A and come to Gaussian integrals that can be 
treated without difficulties. For the problem considered the propagator K is proportional to the average of the square 
of the field (|A/c| ). In terms of the functional integral this quantity is written as 



/ I A|' exp(-/3Fe//[A, A*])DADA* 
J exp{-f3Feff[A,A*])DADA* 



(3.7) 



here (3 = 1/T and Feff[A, A*] is the effective free energy functional. We have chosen the parameters in such a way 
that the grains are zero-dimensional. Therefor e, it is sufficient to integrate over the zero space harmonics only, which 
means that the field A in the integral in Eq. (3.7) does not depend on coordinates. In the quadratic approximation 
in the field A the free energy functional -Fe// includes two different terms 






p(2) 



(3.8) 



where F^J^ describes the superconducting fluctuations in an isolated grain and F^rK takes into account tunneling from 
grain to grain. For the first term we obtain after standard manipulations 



i^S[A,A*]=F5]ll/|A|-r Y. 2C(zs„)||A(f7,)r 

Ofc \ 2e„>|Ofc| 



(3.9) 



where the function C(i£„) is defined in Eq.(3.6) and V is the volume of a single grain. In the limit of low temperatures 
T <C £{){H) the sum over the frequencies e„ in Eq. (3.9) can be replaced by the integral and we reduce the functional 

FJ]^ [A, A*] to the form 



fW[a,a* 



^i:(.» 



Ofc 



£o{H) 



\n. 



Ao 



|A(f]fe) 



(3.10) 



Close to the critical magnetic field He destroying the superconducting gap in a single grai n the energy of the first 
harmonics £o{H) is equal to the BSC gap at zero temperature Aq. This means that Eq. ( ^.1C| ) can be w ritten in 
this case in the region T <^ Aq. Near He small frequencies Vlk are most important and one can expand Eq. ( 3.10 ) in 
powers of the small parameter ilk/^o- Then, we obtain 



^S[a,ai4e^ 



go(g) 
Ao 



\Jh\ 
Ao 



iA(f]fc)r 



(3.11) 



At strong magnetic fields H ^ He, one should use the more general formula, Eq. (3.10) 
eff 



The term -^^V^ [A, A*] describing the tunnchng includes three different contributions represented in Fig. 5 
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FIG. 5. Diagrams describing FKr. 
The analytical expression F^tt' [A, A*] corresponding to the first diagram Fig. 5 can be written as 

F^f/ [A, A*] = -ty{27:i.orfJ2^A,A*+c.c.)V' / ^!^l^rE G(z£„pOG(-z£„pOG(ze„p,)G(-z£„p,)C2(z£„) 

J {2n) ^^ 

(3.12) 



p(21) 



^,3 



Writing Eq. (B.12) we put Clk — in the expression for the Cooperon G and in the Green functions. This is justified 



because the energy F^A [A, A*] is already small because it includes the parameter J {5 / Ea{H)) that is assumed to be 



small, where J is the dimensionless conductance of the system specified by Eq. ( ^.11 ). Next terms of the expansion 
are of the order of J{S/£o{H)){ilk/£o{H)) and can be neglected for small ilk- The second and third diagrams in Fig. 



5 are equal to each other and have the opposite sign with respect to the first diagram. For simpHcity we assume that 
the granules are packed into a cubic lattice. Using the momentum representation with respect to the coordinates of 
the grains and taking into account all diagrams in Fig. 5 we reduce the free energy functional F^J. [A, A*] to the 
form 



p(2) 



E 

i=l 



F^fj [A, A*] = (8/3^)(l/<5) ^ J{S/£oiH)) (1 - cosq^d) |A| 



where q is the quasi- momentum and d = 2R. Eq. ( 3.13| ) is written in the limit 

J-^£o{H)/5 
The inequality ( [3.14 ) is compatible with the inequality ( 2.12| ) provided the inequality 

£o{H) » S 



(3.13) 



(3.14) 



(3.15) 



is fulfilled. If £q{H) ^ Aq, the condition, Eq. ( 3.15 ), is at th e same time the condition for the existence of the 
superconducting gap in the single granule. The inequality ( t3.14[) allows us also to neglect influence of the tunneling 
on the form of the Cooperon, so we use for calculations Eq. ( |3.6D . 

(2) 

Writing the previous equations for F^A [A, A*] we neglected the influence of magnetic field on the phase of the 
order parameter A. In other words we omitted the phase factor exp(z- J A{r)dr). The effect of the magnetic field on 
the phase will be discussed in details later in Sec. Vl. 

Although the final result for the correction to the DOS can be written for arbitrary te mperatures T and magnetic 
fields H, let us concentrate on the most interesting case T ^ Tc, H > He- Using Eqs. ( p/7 3.1C , 3.13 ) we obtain for 
the propagator of the superconducting fluctuations K (iflk.q) 



K{ink,q) 



In I -^^ ) + V{(i. 



77(q) = (8/3^) J2 J iS/^o{H)) (1 - cos q,d) 



(3.16) 



j=i 



The pole of the propagator K{ink, q) at q — 0, flk — determines the field He, at which the BCS gap disappears in 
a single grain. From the form of Eq. (3.16) we find 



£o(ffc) = Ao. (3.17) 

Th e result for He, Eqs. (3.5, 3.17), agrees with the one obtained long ago by another mcthodEJ. We can see from 
Eqs. ( 3.16|, |3.17J) th at the term 77 (q) describing tunneling is very important if H is close to He- 

Eqs. ( |3.6| ), ( |3.16 ) give the explicit formulae for the functions C and K and allow us to calculate the correction Sv 
to the DOS . The analytical expression for the diagram. Fig. 2a, reads as follows 

Siy{ien) = {l/7T)l/{2TTi^oTf [ ^-^T Y^ K {iQ k , €[) [^C^{ien,ink~ien)G^ien,p)G{tnk-ien,p) (3.18) 

J (27r) ^^ J 2tt 



Eq. ( 3.18| ) contains integration over the momentum in the single grain p and the quasi-momentum q. First, we 
integrate over the momentum p and reduce Eq. (3.18) for £„ > to the form 



Sv{ien) = y^ / K{ink,q)C'^ {ien,i^k 



iC.f-1 



Ofc<e 



Sq_ 



(3.19) 

e. At low 



After calculation of the sum over O^ in Eq. (3.19), one should make the analytical continuation iSn 
temperatures, it is suf ficient to find the correction to the DOS at zero energy Sv = dv (0). 

Remarkably, Eqs. (3.16-3.19) do not contain explicitly the mean free time t. This is a coiisequence of the zero- 
harmonics approximation, which is equivalent to using the random matrix theory (RMT) Ref.tJ. (The parameter r 
enters only Eq. (|3.5|) giving the standard combination Sq (H) describing in RMT the crossover from the orthogonal 



to the unitary ensemble). This justifies the claim that the results can be used also for clean grains with a shape 
providing a ch aoti c electron motion. 

Using Eqs. ( p^ , |3. 16 -3. 19) one can easily obtain a n exp licit expression for odos for H — He <^ He- In this limit, 
one expands the logarithm in the denominator of Eq. ( |3.16|) and neglects the dependence of C on g„ a nd ilk because 
the main contribution in the sum over flk comes from i7fc '^ £o{H) — £q{Hc) ^ Aq. Using Eq. (2.14) the result for 
S<Jdos = o'o ^ o-DOS can be written as 



d<7 



DOS 



CTO 



2S 



-TT-i <ln^(q) >,, r/Ao<^ 

^<77-i(q)>„ r)«r/Ao«l 



(3.20) 



2ir/d 



7?(q) = ^(q) + 2/i, < 



>q^V 



{...)dq/{27rf 



where h — {H — He) /H^ and V^ is a volume of the single grain. We see that the correction to the conductivity is 
negative and its absolute value decreases when the magnetic field increases. The correction reaches its maximum at 
H —>■ He- At zero temperature and close to the critical field H^ such that J((5/Ao) S> h, the maximum value of 
So'DOs/ci'o from Eq. ( 3.20| ) is 
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(3.21) 



In the limit J{S/Aq) <$:. h < 1, one can expand the logarithm in Eq. ( 3.20| ). Then, taking h ^ 1 the correction to the 
conductivity at zero temperature can be estimated as 



5(tdos 



CTO 



(3.22) 



Schematically, the suppression of the DOS due to the superconducting fluctuations is shown in Fig. [6[ 

V/Vo, 



/Ao 



FIG. 6. Suppression of DOS due to superconducting fluctuations. 



As temperature grows, the correction to the conductivity due to the reduction of the DOS can become larger and 
reach for T ^ Aq and J((5/Ao) 3> h the order of magnitude of J~^ . 
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(3.23) 



In the limit h ^ I ^ J(6/Ao) and at temperature T ^ Aq this correction can be estimated as 






_S_ 
A'o 



(3.24) 



We see from Eqs. ( p.21 -3.24) that the corrections to the conductivity are smaller than unity provided we work 
in the regime of a good metal, Eqs. (2.5, 2.12), so the diagrammatic expansion we use is justified. Indeed, we can 
neglect the corrections of higher orders. For example, the diagram shown in Fig. 7 has the additional small factor of 
(S/Aq) Inf] at T/Aq < 77 and S/{Ao^) at r) < T/Aq. 




FIG. 7. An example of high order corrections to DOS. 



The correction to the conductivity calculated in this Section could become comparable with ctq when J ~ 1. 
However, such values of J mean that we would be in this case not far from the metal-insul ator transition. Then, we 
would have to take into account all localization effects. Epr values of J ~ 1 one can use Eq. ( 3.20| ) for rough estimates 
only. Apparently, the parameters of the samples of Ref.El correspond to the region J ^ 1, (5/Ao ^ 1/3. 

In the limit of strong magnetic fields H :^ He the correction to uo can still be noticeable. In this case we can 
use Eq. (3.1C) as before but, with a logarithmic accuracy, we can neglect the dependence of the superconducting 
propagator, Eq.(3.16), on Qk and on the tunneling term. Then we obtain finally 



Sa 



DOS 



/ao = -(1/3) {S/£o{H)) In"! {£o{H)/Ao) 



(3.25) 



Eq. (3.25) shows that in the region H ^ He the correction to the conductivity decays essentially as daoos ^ H^^ . 

Let us emphasize that the correction to the conductivity coming from the DOS remains finite in the limit T ^ 0, 
thus indicating the existence of the virtual Cooper pairs even at T = 0. 

In the region of not very small h 3> JS/Aq, we can neglect the tunneling term F^A [A, A*] in the free energy 
functional F^ff [A, A*]. Then, we can write the correction to the conductivity in a rather general form. The super- 
conducting propagator K can be written in this case as 



K{ink) = -i^o^ln 



£o{H) + pK 
Ao 



(3.26) 



Using Eq. ( 3.19 ) for the correction to the DOS and Eq. ( 2.14 ) we obtain for the correction to the conductivity at 
zero temperature 



Sapos _ _ 1 / S 
(Jo 3 IAq 



X 3 \Ao' 



(3.27) 



In the limit ft <C 1, we reproduce with logarithmic accuracy Eq. (3.21), whereas in the opposite limiting case /i 3> 1 
we come to Eq. ( 3.25| ). 



In order to calculate the entire conductivity, Eq. (2 



), we must investigate the AL and MT contributions (Figs. 2c 
and 2b). In conventional superconductors near Tc, these contributions are most important leading to an increase of 
the conductivity. In the granular materials, the situation is much more interesting. It turns out that both the AL and 
MT contributions vanish in the limit T ^ at all iJ > i/c and thus, the correction to the conductivity come s from 
the DOS only. So, at low temperatures, estimating the total correction to the classical conductivity cto, Eq. ( 2.10 ), 
one can use the formulae of this Section. 



IV. ASLAMASOV-LARKIN CORRECTION TO THE CONDUCTIVITY 

The Aslamasov-Larkin (AL) correction to the conductivity ctal originates from the ability of virtual Cooper to 
carry an electrical current. In contrast to the one-electron tunneling determining (TdoSi the probability of tunneling 
of the Cooper pairs from one grain to another is proportional to f*. The quantity aAL is related to the response 



function Q as 



TAL 



Q^^/(-ic^),c^^O 



where the diagram for the Q^^{iujn) is represented in Fig. 2c. Calculating integrals corresponding to this diagram 
we may put Q^ = in the electron loops, because all singularities in the vicinity of the transition point are contained 
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in the propagator of the superconducting fluctuations K^irikp. The analytical expression for the diagram in Fig. 2c 
has the form 



where -Bi(0, q) corresponds to one electron loop. The analytical expression for this loop reads 



i3i(0, q) = — 4i / sinq^dcos(gi — q'^)d 

J f27r 



(27r)^ i2TTiyoTf 



(4.1) 



(4.2) 



d^Pid^P2 
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T^G{ien, pi)G(-i£„, pi)G(ie„, p2)G(-i£„, P2)G^(ie„) 



The functions singed and cos(gi — q'^^d in Eq. (4.2) correspond to the current and tunneling vertices, respectively. Eq. 
(4.1) is obtained by considering four different types of AL diagrams that are obtained from each other by permutations 
of the current and tunneling vertices. Summing over the spin of the electrons we get the additional factor 2. As in the 
preceding Section we calculate the impurity vertices neglecting the tunneling term, which is justified if the inequality 
( |3.14 ) is fulfilled. Integrating over the momenta pi and P2 in Eq. (4.2) we reduce the functions i?i(0,q) to the 
following form 



^l(0,q) = 2 J 



d3q' 



-^■J c ^u\ 4' / sinq,dcos(qi - q^)d—— 
TT^ toiH) J (27r 



(4.3) 



To calculate the response function Q for real frequencies uj one has to make an analytical continuation from the 
Matsubara frequencies w„. This can be done rewriting the sum over Qk iihEq- (4.1) in a form of a contour integral 
that allows to make the continuation ia;„ -^ ui + iO. As a result, we obtainQ 



tJ2k (i^k, q) K {iflk - lUJn, q) 



(4.4) 
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+ 00 /■ + 00 

coth ;^(i^^(£) - K^{e))K^{e - Lj)de + / coth ^(^^(e) - K^{e))K^{e + Lo)de 



where K^{K ) is retarded (advanced) superconducting fluctuation propagator. Expanding Eq. (4.4) in lo we keep 
the term that remains finite in the limit w — > and the linear one. The zero order term cancels with the contribution 
of a diagram schematically represented in Fig. 2c but containin g in stead of the current vertices the tunneling ones. 
The latter contribution originates from the second term in Eq. (^J). The linear term giving the dc conductivity can 
be written as 



tY,K (irifc, q) K {iQk - luJn, q) 



2u£i{H) 



coth Trde 

2T (^2 + £2)3 



(4.5) 



where 77(q) is defined in Eq. (3.16). Using Eqs. (4.1, 4.3,4.5) the fluctuational contribution gal to the conductivity 
is reduced in the limit T <^ Tr. H — Hr -^ Hr to the form 



— = Y-^^ E < ^ (q) [^'(q)]-' s^\d >, 



CTO 



(4.6) 



^0 ,-_i 



where A (q) = AttT (3Ao?7(q)) ^ for T < Agry and A (q) = 1 for AqV <$:T <Tc. 

Using Eq. ( |4.6| ) we can estimate the quantity (tal/<^o- 

Let us consider fir st t he limit T <C Aq?]. In this region, provided the inequality h <C JS/Aq is fulfi lled, the main 
contribution in Eq. (4^) comes from small q. Calculating the integral over q we obtain from Eq. (4.6) 



aAL 32n^_ jS^T^ 



O-Q 
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(4.7) 



From Eq. (4.7) we can see that at low temperatures the AL correction to the conductivity is proportional to the 
square of the temperature ctal ^ T^ and vanishes in the limit T —f 0. 
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Let us co mpar e the AL co rrection with correction due to suppression of the DOS considered in the previous Section. 
Using Eqs. (4J) and ( 3.21 ) we obtain 
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(4.8) 



We see from Eq. (4.?) that at T ^ Ao?y the AL correction is small \(Jal/5'^dos\ ^ 1- This means that the AL 
contribution cannot change the monotonous increase of the resistivity of granulated superconductors when decreasing 
the magn etic field. At very strong magnetic fields H ^ He, we can neglect the second term in the denominator of 
Eq. (4.7). Then the AL correction can be estimated as 



<^AL 
CTO 



<52t2 
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(4.9) 



To compare this result with the correction due to the DOS we should use Eq. (3.25) that was also derived at strong 
magnetic field H :^ He- 
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Using Eq. (3^) for £q{H) we can see from Eq. ( 4.10 ) that \(JAL/S(JDOs\H-»Ha ^ H^^. 

N ow let us consider the region of temperatures not far from the critical temperature Aq?) <C T < Tp. Using Eq. 
(4.6) we have 
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(4.11) 



In the limit J((5/Ao) 3> h, the main contribution in Eq. (4.11) comes from small q and the contribution aAL can be 
written as 



CTO 
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(4.12) 



We see that the AL contribution grows when approaching the critical field Hr- In order to calculate the AL 
correction we used the perturbation theory. Therefore, the region of the validity of the results obtained is described 
by the inequalities aAh/fyo <C 1 or /i :^ J'^{Aq/5). 

Now let us compare the AL correction with Sa^os- From Eqs. (4.11, 3.23|) we obtain 
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(4.13) 



Eq. (4.13) is correct only at fields close to the critical field, such that the inequality h ^ J~^{Ao/6) is fulfilled. In 
this region the total correction to the conductivity is positive, which means the resistivity decays, when approaching 
the critical magnetic field He- In the case of a strong magnetic field H ^ He, we can neglect the second term in the 
denominator of Eq. (4.11) and then we obtain 



O'AL 

(To 



0(§ 



(4.14) 



that is aAL ~ H^^ at iJ > i/c and T -- Aq. 

To understand the behavior of the total conductivity of the granulated superconductors in this region we should 
consider also the Maki-Thompson contribution and this will be done in the next Section. 

V. MAKI-THOMPSON CORRECTION TO THE CONDUCTIVITY 

Another contribution usually increasing the conductivity is the Maki-Thompson (MT) contribution represented 
in Fig. 2b. Again, we can put fife = in the electron loop, because characteristic frequencies in superconducting 
fluctuation propagator K are of the order fi^ ^ A077 <C So- The analytical expression for this diagram reads 



12 



3 /. 73 

Q^'^itUn) = 2/3 V / — ^T Vx(zr!fc,q)i32(*t^„,q) 
. , J {2n) „ 



where B2{ioJm q) is a function describing the contribution of the loop. This function can be written as follows 



(5.1) 



B2(iw„,q) = 4 / sing-dsin(gj - q[)d- — 
J (2ti 



{2i,f {2^v^rf 



d^Pid^P2 
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r^G(-ze„,pi)G(ie„,pi)G(-ie„ + icj„,pi)G(ie„ - icj„,p2)G(ie„)G(-i£„ + icj„,ze„ - za;„) (5.2) 



where, as before, pi and p2 stand for the momenta in the granules and q, q' are quasi-momenta. Eq. ( |5.lD includes 
the contribution of two diffe re nt M T diagrams and summation over spins. In order to make the analytical continuation 
iwn — > o; + iO in Eqs. (5.1), ( |5.2| ) it is convenient to rewrite the sum over e„ in the form of the following contour 
integral 

T^ ... = -^( /c^ tanh ^G^("^, Pi)G^(^, Pi)G^(-z + ic^„, p2)G^(z - ^c^„, p2)G^(z)G«(z - iujn)dz + 



/p^ tanh ^G^(-z, pi)G^(z, pi)G^(-0 + iu;„, P2)G'^(2; - iw„, p2)G^(z)G^(z - icj„)dz + 

/c, tanh :^G«(-z, p^G^^lz, pi)G«(-z + il^„, P2)G'^(z - ic^„, p2)G^(z)G'4(z - zu;„)dz), (5.3) 



where the contours Gi, G2, G3 are shown in Fig. 8 



" ^1 



C. 



c. 



FIG. 8. Contour of integration. 
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C2 
C3 



As usual, the MT diagrams have both regular (contours Gi, G3) and anomalous (contour G2) part. Near T^ and at 
very low magnetic field the anomalous part can be very large. In the limit iJ — ^ £L, it can even diverge and become 
larger that the AL correction giving a positive contribution to the conductivityt2l. However, in the limit of high 
magnetic fields and T ^ 0, the situation is less intriguing. It turns out that for the problem considered, the absolute 
values of the regular and anomalous parts are equal in this limit but these contributions have the opposite signs. 
Making the analytical continuation in Eq. (5.3) and integrating over the momenta we obtain in the lowest order in lo 
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(5.4) 



where the energy Eq = ^o (H) is given by Eq. (3^). From Eq. (5^), we see that at T = the MT contribution 
vanishes, which corresponds to the cancelation of the regular and anomalous parts. At low but finite temperatures 
T, the final result for the MT contribution can be written as 
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(5.5) 



where B (q) = -tt^i ln?7(q) for T < Aq?? and B (q) = 2T (Ao?7(q)) ^ for Aq?? <.T <Tc 

Let us estimate the MT contribution in the different limiting cases. At low temperatures T ^ Apry and h ^ J(5/Ao 
we have from Eq. (5.5) 
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(5.6) 
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Comparing Eq. (^.6|) with Eq. (3.21) we come to the following estimate 
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If the magnetic field H is not very close to He, such that JS /Aq ^ /i < 1, we obtain from Eq. (5.5) 
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(5.8) 



We can see from Eqs. (5.7, p^) that the MT contribution (Jmt is proportional at low temperatures to T^, which is 



the same temperature dependence as that for the AL contribution. This means that, at sufficiently low temperatures, 
the MT contribution is sm all, <j mt I Sapps <^ 1. Thus, we conclude that, in this region, the main correction to the 
classical conductivity, Eq. (2.10) comes from the correction to the DOS, Eq. (3.21). The latter correction is negative, 
so the resistivity of the granulated superconductors exceeds its classical value. 

It is interesting to compare the AL and MT corrections. In the limit T ^ 0, and h ^ JS/Aq, we obtain using 
Eqs. 1^ and (|^) 
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Eq. (5.9) shows that in this region. The AL cont ribu tion is larger the MT one. Let us consider another case of not 
very low temperatures, Aqt? <C T < Tc- From Eq. (5.5) we have 
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(5.10) 



Eq. ( 5.10 ) gives the possibility to compare the MT contribution with the DOS in this temperature interval. Recalling 
Eq. (3.23) we obtain 
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Eq. (5.11) shows that at not very low temperatures the MT contribution has the same order of magnitude as the 
contribution due to the reduction of the DOS. At the same time, the AL contribution in the region Aqt) <^ T ^ Tc 
can be consid erably large r tha n both the MT and DOS contributions. 
From Eqs. ( [4.12| ) and ( |5.10| ) we can see that at T < T^ 
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Eqs. (5.11,5.12) show that, at not very low temperatures, and not far from the critical field the AL correction to the 
conductivity is the most important. This means that approaching the transition in this region the resistivity decreases, 
which is in contrast to the behavior at very low temperature where the correction to the resistivity is determined 
entirely by the contribution to the DOS and is positive. 

To conclude the last two Sections we emphasize once more that the temperature and magnetic field dependence 
of UAL and umt is rather complicated but they are definitely positive. The competition between these corrections 
and (Tdos determines the sign of the magnetoresistance. We see from Eqs. (4.6, 5.5) that both the AL and MT 
contributions are proportional at low temperatures to T^. Therefore the (tdos in this limit is larger and the magne- 
toresistance is negative for all He- In contrast, a,t T ^ Tc and close to He, the AL and MT corrections can become 
than <7dos resulting in a positive magnetoresistance in this region. Far from He the magnetoresistance is negative 
again. 



VI. THE CRITICAL FIELD Hc2 IN THE GRANULATED SUPERCONDUCTORS 

In the previous Sections we considered transport in granulated superconductors at magnetic fields H not far from 
the field He- The field He is the field destroying the superconducting gap in a single isolated granule. We have seen 
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that the m ain contribution due to the superconducting fluctuations comes from the correction to the density of states, 
Eq. ( 3.20 ), and this correction remains finite in the hmit H ^ He- But is the field He a critical field in the system of 
the granules coupled to each other by tunneling? If it were a critical field what would happen aX H < Hcl Would the 
system be macroscopically the superconductor or normal metal? Or, may be, this wo uld be a new st ate of matter? 

To answer these questions we should derive the effective action F^ff [A, A*], Eqs. (3.8, 3.10, 3.13), more carefully 
than it has been done in Sec. [II. Namely, until now we considered only the effect of the magnetic field on the electron 
motion inside the grains, neglecting its influence on the correlation of the phases of the order parameter Ai of different 
grains. However, to understand whether the system is macroscopically superconductor or not, we must consider the 
macroscopic motion and thus, the effect of the magnetic field on the phase correlation. 

In this Section we come back to the derivation of effective action fe// [A, A*] taking into account the influence 
of magnetic field on the phases of the order parameter. First, we consider this problem qualitatively and then, 
quantitatively. It is clear that if the magnetic field is strong enough, it induces macroscopic currents that finally, at 
a field iJc2 j destroy the superconductivity. 

Let us estimate the critical magnetic field Hc2 using the Ginzburg-Landau theory. We assume that the granulated 
system under consideration is in the macroscopically superconducting state. Then, the Joscphson part of the free 
energy in the coordinate representation, not too close to the Hc2 , can be written in the form 
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E^..^ 



A,)^ 



A^ 






(6.1) 



where Ej = JAq is the Josephson energy, J ^ 1 is the dimensionless conductance, Eq. ( 2.11 ), Ao is the BCS gap 
at zero magnetic field and i? is a radius of the single grain. The gradient expansion in Eq. (6.1) was done under the 
assumption that the magnetic flux through one grain is smaller than the flux quantum (pQ. 

In the conventional Ginzburg-Landau free energy the co effic ient in front of the gradient term is proportional to 
the square of the coherence length ^. Therefore, using Eq. (3.1) we can extract the macroscopic coherence length ^. 
Recalling that the free energy of a single grain is i/VA'^ , where 1^ is a volume of single grain, we obtain 



e - JR\5/A^) 



(6.2) 



If the conductance J is large enough the behavior of the granulated superconductors is the same as in a bulk sample 
with the effective coherence length ^. Now we estimate the critical magnetic field Hc2 th at de stroys the superconduc- 



tivity in the system as Hc^^^ « (po- This field is different from the field He- Using Eq. ( 3.17| ) we can compare these 
two fields. The result for the ratio of these two fields is 




(6.3) 



Eq. (6^) for Hc^ is valid provided Hc^ ^ -ffc, which corresponds to the inequality J » {Aq / S) y^^o/R- However, we 
are interested in the opposite case when 



1« J«(Ao/<5)v/W^ 



(6.4) 



This contradict to the assumption made and means that, in the region specified by Eq. (6.4), the critical field Hc2 is 
close to the field He and \Hc — Hc^] /He can be considered as a small parameter. 

Now let us calculate the critical magnetic field H^^ more rigorously taking into account the influence of the magnetic 
field on the macroscopic motion. The magnetic field results in an additional phase factor exp(i- J A(r)dr) in the 
s upercond ucting field A^ in Eq. ( |3.12 ). We assume that this magnetic field is not far from He determined by Eqs. 
( |3.5| , |3.17|) . The field He is of order (po/ RS,, which means that, in the limit under consideration i? ^ ^, the magnetic 
flux through the elementary cell of the lattice of the granules is small. Therefore, w e can expand the functions 
exp(z| / A(r)dr) in A and write gradients instead of the finite differences of A^ in Eq. ( 3.12 ). Essential frequencies 
ilfc are also small and the free energy functional in such a continuum approximation takes the form 



F[A,A*] = \Y. 
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(6.5) 



The critical magnetic field iJc2 c an be found writing the propagator of the superconducting fluctuation K corre- 
sponding to the free energy, Eq. (3.5). Making Fourier transformation of the function A in the eigenfunctions of the 
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operator entering Eq. (6.5) and calculating Gaussian integrals we obtain for the propagator Jf„ (0, qz) in the spectral 
representation at Clk — 



Kn{0,qz 
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(6.6) 



where qz is a component of the quasi- momentum parallel to the magnetic field and n = 0, 1... are the number of the 
Landau levels. To calculate the critical magnetic field Hc2 we should consider poles of the superconducting propagator. 
Taking the lowest Landau number n = and putting qz ~ obtain the following equation determining the critical 
field Hr. 
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Expanding the first term in Eq. ( p.7| ) near He we find the critical field He 
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(6.8) 
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Eq. (|6.8D shows us that the critical field Hc2 is close to the field He s o long as {JS /Aq) {R/£,) <S^ 1, where ^ ^ {£,ol) 
is the coherence length in the superconducting grains. Eq. (3.f4) and the assumption that the grains are zero- 
dimensional guarantee the fulfillment of this inequality. For the ballistic motion of electrons inside the grain (the 
radius of the grain R is of the order of the mean free path, i? ~ Z) the critical field Hc2 can be written in the form 



He- He, 40 f 5 \ R 

He ^2 yAoJy^o 



(6.9) 



Eq. (6.8) is the main result of this Section. Below the magnetic field He2, one should add in Eq. (6.5) a term quartic 
in A, which gives a non-zero order parameter A. Thus, at field H < Hc2 the granular system is in the superconducting 
state. 

In order to understand the behavior of the resistivity as a function of the magnetic field H in the region Hc2 < 

H < He v/e can consider the quantity -Mj ( ^)- We can use as before Eq. (B.f£) but now, calculating the integral, 

we should make the following replacement 
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where qz is a component of the quasi-momentum parallel to the magnetic field. 
The correction to the DOS at T = takes the form 
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The main contribution to the correction 5v{Q) in Eq. (3.1C) comes from the term with n = Q. Using the fact that 



5a / aa — 25v /v^^ taking the first derivative with respect to the magnetic field and finally, integrating over the frequency 
il and quasi-momentum qz we obtain 
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(6.11) 



where He2 is given by Eq. (6.8). We can see from Eq. (|6.11) that the value -Mj- {^) diverges when H approaches 



He2 ■ Thus, the critical field i?c2 is characterized by the infinite slope on the dependence of the resistivity on the 
magnetic field. This property might help to identify this field on experimental curves. 
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VII. ZEEMAN SPLITTING. 

In our previous consideration we neglected interaction between the magnetic field and spins of the electrons. This 
approximation is justified if the size of the grains is not very small. Then, the critical field H°'^ destroying the 
superconducting gap is smaller than the paramagnetic limit g^H^ = Aq and the orbital mechanism dominates the 
magnetic field effect on the superconductivity. However, the Zeeman splitting leading to the destruction of the 
superconducting pairs can become important if one further decreases the s ize of the grains. 

Let us discuss now the effect of Zeeman splitting. We can rewrite Eq. (1.4) for the ratio of orbital H"^ magnetic 
field to the Zeeman magnetic field in the following form 



H^ 
H? 



1/3 



(7.1) 



where Sc ~ l/{i'R'^), Re ~ C(PoO ^ ^-^^d ^ — VCo^- To understand whether the Zeeman splitting is important for an 
experiment we can estimate the ratio (5c /Aq and compare it with the proper experimental result. We find easily 



k = ^'"lTo 



(7.2) 



which shows (5e/A o ca n be both smaller tha n {S /Ao)exver. and larger depending on the values of I and ^q. Using the 
result for Sc, Eq. ([7.2|), we can rewrite Eq. (7.1) as 



AT 



_S_R 



1/2 



(7.3) 



In the experimental, both the mechanisms are in principle important. One can come to this conclusion using the 
fact that the Zeeman critical magnetic field is H^ « 3.5T and this is not far from the peak in the resistivity at the 
field H « 2.5T. Below, we consider the corrections to the DOS and conductivity due to the Zeeman mechanism. We 
will see that at temperature T = these corrections can be of the same order of magnitude as the correction due to 
orbital mechanism. 

Let us calculate first the critical magnetic field Ho destroying the superconducting gap in a single grain taking into 
account both the orbital and Zeeman mechanisms of the destruction. The Green function for the non-interacting 
electrons in this case is 



Gf^{ien,p) 



1 



i£„ - e(p) ± Ez/2 + (i/2T)sgn{en) ' 



(7.4) 



where Ez = gfJ^sH is the Zeeman energy. Including the interaction between the magnetic field and electron spins we 
obtain the following form of the Cooperon 



C(ie„, -i£„) = 27riyo(|2e„| - iEzsgn (e„) + £a{H)y 



(7.5) 



Repeating the calculations of Sec. [II with the modified Cooperon, Eq. (L5), we find at T ^ Tc the new critical 
magnetic field Hq 



S^Ho) + EliHo) = Al 



(7.6) 



In the limit of a very week Zeeman splitting, when the o rbital mechanism is more important, Eq. (7^) reproduces the 
previous result for the critical magnetic field Eq. (3.17). In the opposite limiting case, when the Zeeman mechanism 
plays the major role, we obtain Ez = Aq. In g ener al case, when the both mechanisms of the destruction of the 
conductivity are important, one should solve Eq. (7.6) and the result reads 



^0 — He 




1 fgfisH, 



1 fgfisHc 



1/2 



2 V Ao 



4 V Ao 



(7.7) 



Using Eq. (3.17) for a critical field He we can estimate the ratio fisHe/Ao = {S,o/R)iPoR) (pol) ^1- If this 



parameter is small (the orbital mechanism is more important) the critical field Hq is close to the field He 
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Ho^Hjl-'y^^'^''^ 



4 V Ao 



(7.8) 



When Zeeman splitting is more important then orbital mechanism then from Eq. (7/7) we obtain g^sHo — Aq. This 
is the point of the absolute instability of the paramagnetic state. At gusH < Aq, the superconducting state is the 
only stable one. 

Let us calculate the correction to the DOS taking into acc ount both the Zeeman splitting and orbital mechanism 
of the suppression of superconductivity. Using Eqs. (7.4, 7.5) we obtain for the superconducting propagator K 



i^(^^.) = -crM]in(i^^^^L, 



nk\y 



El 



'7(q) 



(7.9) 



In the region h = {H — Ho)/Ho <^ 1, where iJg is given Eq. (7.7), expanding the logarithm in the superconducting 
propagator we obtain 



Kiiflk) = -Vq' 2h 



Ao 



7?(q) 



(7.10) 



Using Eqs. ( 3.19| , 2.14) for the c orrection to the conductivity we obtain the same result as before, Eq. ( ^.20 ), but 
with the new h^. Eqs. ( [j. 21 - 3. 24 ) are correct for in general case. 

In the limit of a strong magneti c fie ld h ^ I we can neg l ect w ith logarithmic accuracy the ri-dependence of 
superconducting propagator in Eq. (7.9). Using Eqs. ( 3.19| , 7^3, ^.14 ) we obtain for the correction to the conductivity 
at strong magnetic fields 



Saoos 



S 



3 \EziH) 



arctan 



f Ez{H) \ 
\ So{H) ) 



\n 



£§{H)+EliH) 



^0 



(7.11) 



From Eq. (7.11) we can see that if orbital mechanism is more important than the Zeeman one, that is if Ez /Sp <C 1, 
we reproduce the previous result for the correction to the conductivity at strong magnetic field, Eq. (3.25). If Zeeman 
mechanism is more important Ez/Sq 3> 1, then we obtain for the correction to the conductivity 



Sa 



DOS 



O-Q 



6 \Ez(H) 



In 



EzjH) 

Ao 



(7.12) 



Eq. ( 4.13 ) has the same structure as Eq. ( p.25| ) but the function 8q{H) is replaced by Ez{H). This changes the 
asym ptoti c behavior at strong magnetic fields because Ez [H] ^ H in contrast to Eq (H) ~ H^. So, we conclude from 
Eq. (|7.12|) that SaDos/<ya ~ H-\ 



VIII. DIAMAGNETIC SUSCEPTIBILITY OF GRANULAR SUPERCONDUCTORS 



In previous sections we have demonstrated that the resistivity of the granulated superconductors grows when 
approaching the superconducting state from the region of very strong magnetic fields. Resistivity is a quantity 
studied experimentally most often. Another quantity accessible experimentally is the magnetic susceptibility. Can 
one observe anything unusual measuring the dependence of the susceptibility as a function of the magnetic field? 

The diamagnetic susceptibility of a bulk sample above the critical temperature Tc in a weak magnetic field has been 
studied long agots. In this Section, we want to present results for the diamagnetic susceptibility x of the granular 
superconductors in the opposite limit of strong magnetic fields and low temperatures. As will be shown below, the 
fiuctuations of virtual Cooper pairs always increase the absolute value of the diamagnetic susceptibility x of the 
granulated system. In contrast to the conductivity, the diamagnetic susceptibility is determined mainly by currents 
inside the granules and is finite even in isolated granules. Therefore, the fact that the virtual Cooper pairs cannot 
move from grain to grain, which is crucial for the conductivity, is not very important for the magnetic susceptibility 
and the latter does not show a non-monotonic behavior characteristic for the resistivity. 

To derive explicit formulae, let us consider first the limit of very low temperature T <^Tc and strong magnetic field 
H — Hc <C He- The effective free energy functional F [A, A*, H] in the quadrati c approximation in the order parameter 
A has been already obtained for this case and is given by Eqs. (3.8, 3.11, 
be calculated using the standard relations 



3.13). The diamagnetic susceptibility x can 
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X^-{l/V){d^F/dH^), 



(8.1) 



where the free energy F has the form 

F = -Thi ( I dAdA*exp[-(3Feff{A,A*,H)] 



(8.2) 



where, as before, V is the volume of a s ing l e gra i n and the averaging (...)g is specified in Eq. (3.20). Calculating the 
derivative in Eq. (8.1) and using Eqs. (3.8, 3.11, 3.13, 8.2) we find for the diamagnetic susceptibility 



}'1^{''-''' 



He I ^k 

2 An 



37r 2 



— j|](l-cosg.d)j 



(8.3) 



At very low temperature T/Aq <^ h = {H — He) /He, we can replace the sum over frequency flk by an integral. As a 
result we obtain 



■5^H^ 



5^-^(^A^jE(i--«'?^^) 



(8.4) 



where xl = e^uo/127r^?ic^ is the Landau diamagnetic susceptibility. In the limit h 3> J{6/Ao) tunneling between 
grains is not important and the diamagnetic susceptibility of the granulated system is equal to the susceptibility of a 
single grain. If the motion of electrons inside a grain is more or less ballistic (/ w R) the result takes the form: 



xiH)--^XL(j^ 



-10- 



He 



H -He 



(8.5) 



From Eq. (8.5), we can see that the fluctuation- induced diamagnetic susceptibility x can appreciably exceed the value 



XL due to the Landau diamagnctism. In order to probe the diamagnetism due to the superconducting fluctuations 

experimentally one can measure the field-dependent part of susceptibility and compare it with Eq. (8.5). 

In the limit h <C T/Aq, it is sufficient to take into account only one term with flk = in the sum in Eq. ( ^.3| ). 
Then, the result for the diamagnetic susceptibility x can be written in the form: 



X- 



I 
5"^^ V R 



T 

A^ 



37r V Ao 



^(1 - cos Qid) 



(8.6) 



j=i 



In the limit when h ^ J{S/Ao) and I w R ,Eq. (8.6) can be simplified and one comes to the following expression 

3"0' 



X{H) 



n- T 
:XL- 



He 



Ao \H-He 



7) 



Eqs. (8.6, 8.7) show that the diamagnetic susceptibility x diverges in a power law when magnetic field H is close to 
the field He but the powers are different for the two different regi ons of the fields. 

At zero temperature T = 0, as we can see directly from Eq. (8.4), the susceptibility x remains finite even in the limit 
H -^ He- However, as we have discussed in Sec. VI, the field He does not correspond to any phase transition. The 
transition to the superconductivity occurs at a lower field i?c2- So, it is interesting to consider the critical behavior 
near i?c2- 

Proper calculations of the diamagnetic susceptibility in the region He2 < H < He and at te mpe rature T <C Aq can 
be carried out without any difficulty. The free energy functional is given in this case by Eq. ( |6.5[ ) and we obtain for 
the diamagnetic susceptibility x 



— ^E 



He 



dqz 

27r 



H-He,+ 



He l^fcl 

2 Ao 



.. 2^2 JS 2 



In the limit T ^ we replace the summation over flk by integration. Integrating over the frequency and the 
quasi-momentum g^ we obtain 
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We see from Eq. (8.9) that the diamagnetic susceptibility x diverges in a power law when H —t Hc^ but the power 



1/2 is different from those in Eqs. (S.5 



The case of temperatures close to t 



3.7) describing the behavior of the conductivity in the region H > He- 
]he cptical temperature T ~ Tc ^ Tc and weak magnetic fields H -^ He has 
been considered for bulk samples long agollj. So, we present here the result for the diamagnetic susceptibility x of 
the granulated superconductors only in the limit (T — Tc)/Tc and I « R 

X~~XLj-^ (8.10) 



Eq. ( B.IO ) shows that the diamagnetic susceptibility of the granular superconductors near Tc is still larger than the 
magnitude of the Landau diamagnetism. 

In all previous considerations we did not take into account the spin paramagnetism. When the size of the grains is 
large, R^Rc = £,{pol)^^, this effect are small in comparison with the diamagnetism. 

IX. CORRECTION TO CONDUCTIVITY AT |T Tc < Tc 

Effect of superconducting fluctuations on DOS of isotropic bulk samples has been considered in the limit T—Tc ^ Tc 
and H <^ He long aged. The AL and MT contributions to the conductivity were considered in the same limit for 
layered superconductors in RefJl3. Here we want to extend these results to the case of the granulated superconcLuctors. 
Experimentally, only a small increase of the resistivity 5p/ Pq ~ 0.03 — 0.04 has been observed in this regionEl. The 
reason of the reduction of the effect in the vicinity of the critical temperature Tc is that the AL and MT corrections 



are not small in comparison with the correction from the DOS. Repeating the calculations of Sections [II and IV for 
T — Tc^Tc and H = we write the contribution from the DOS and the AL correction as 



ScTDOs 7C(3) S 



ao 27r2 Tc\^^+r,2iq) 



Ta 







(9.1) 



(9.2) 



where 772(9) = 1/3 E JiS/Tc){l - cos q,d). 
i=i 
To calculate the MT correction. Fig. 2b, we should renormalize the impurity vertices taking into account the 
tunneling term. This is necessary, because in the anomalous MT contribution strongly diverges in low dimension if 
the magnetic field is weak. So, the tunnelling from grain to grain can provide convergence of the integrals giving the 
anomalous MT contribution. The Dyson equation for the Coopcron C can be written as 

C = Co + CoSC (9.3) 



where all diagrams for the self-energy S are shown in Fig. 9. The function Co in Eq. (9^) is the Cooperon in a single 
grain. 

^ AW^ ^ , 

-r^— ' : ' 4-Sh-SH-' 

1 J 1 1 J 1 

FIG. 9. Diagrams contributing to the self-energy. 
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Solving Eq. (9.3) we see that the proper propagator for the anomalous part of the MT correction has an additional 
diffusion pole in comparison with the regular part. Therefore, in the limit Tc ^ ??2(9), the anomalous MT contribution 
to the conductivity <J^ is larger than the regular one (J^'^ and their ratio is proportional to Tc/r]2iq)- At the same 

time, the anomalous contribution o■J^^^ is positive, which means that both the AL and MT corrections give a positive 
contribution to the conductivity. Explicit formulae for the regular and anomalous parts of the conductivity can be 
written as 

'^m't _ 7C(3) (5A/ cosg,d \ 
^MT^OOgly/ ^ \ (9.5) 



Eqs. (2_^, 9J., p.2| , j.4 , |9.5| ) describe completely the behavior of the conductivity near T^- We see that the terms 



CTAL and (J^T gi'^iiig positive contributions to the conductivity diverge in the limit T ^ Tc, whereas the terms anos 
and crJv/T reducing the conductivity converge in this limit. Therefore, sufficiently close to Tc, the superconducting 
fluctuations increase the conductivity. A weak magnetic field shifts the critical temperature Tc and one can describe 
also the dependence of the conductivity on the magnetic field. Apparently, far from the transition point one can obtain 
an increase of the resistivity due to the superconducting fluctuations and thus, a peak in the resistivity. However, 
this peak should be small, which correlates with the experimental observation near T^. ft is only the region of low 
temperatures considered in the previous sections where a considerable negative magnetoresistance is possible. 

X. EXPERIMENTS ON AL GRAINS. 

The theoretical study presented in this paper was motivated by the experimental workEI. Let us compare the 
available experimental results with our theory. In the articletJ three samples were studied. We concentrate our 
attention on the samples f and 2, Fig. 4, of that work. 

We analyze the case of very low temperatures T <^ T^ and magnetic fields H > He, where Tc ~ 1.6K is the 
critical temperature for Al grains studied i n the exp eriment and H^ is the critical magnetic field that suppresses the 



superconductivity in a single grain, Eqs. ( |3.5| , p.f7| ). At temperature T ~ 0.3K and magnetic field H ~ AT these 
samples show a large negative magnetoresistance. The resistivity of the sample 2 has the maximum at H = 2.5T and 
the value of this peak is more than twice as large as the resistivity in the normal state (that is, at -ff 3> He, when 
all superconducting fluctuations are completely suppressed). A negative magnetoresistance due to weak localization 
(WL) effects is also not unusual in disordered metals and, to describe the experimental data, its value should be 
estimated as well as the effects of the superconducting fluctuations discussed in the previous chapters. 

The total conductivity of the granular metal under consideration including effects of WL and superconducting 
fluctuations can be written in the form: 

(T = (Jo + Saoos + CTAL + (^MT + SawL (10.1) 

At low temperatures T <C Tc, the contribution Sauos originating from the reduction of DOS due to the formation of 
the virtual Cooper pairs is larger than the contributions ual and umt since the latter vanish in the limit T — > 0. So, 
let us concentrate on estimating the contributions Saoos and SawL ■ 

It is clear that the sample 1 undergoes a metal-insulator transition, which results in the complete suppression of 
the superconductivity. The parameters of the sample 2 that is of the main interest for us are not far from the those 
of the sample 1. Using Eq. (2.II), the value R — 60 A for the radius of the grains, and the value of the resistivity 



po Ki 1.9 X 10~^r2cm we find J « 0.1. 

The small value of J is not in the contradiction with the possibility for the system to be in the metallic phase.|-T|iie 
Anderson metal-insulator transition in granular metals was considered using an effective medium approximationll3ll3. 
The jejjitical point Jc in the present notations is given for the 3D cubic lattice by the equation (Eq. (12.67) of the 
booklia) 
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We can see from Eq. (10.2) that the critical value of Jc = 10~^ is really very small. Therefor, we believe that the 
metal-insulator transition observed in the sample 1 is not a conventional Anderson transition. Apparantly it occtirs 
due to formation of the superconducting gap. Then, the transition can be described following the scenario of Ref.H. 

Why can one be sure that the experimentally the weak localization corrections are small? A similar effect of 
the negative magnetoresistance has been observed in Ref.0 and the authors of that work attributed it to the weak 
localization effects. Could it be that this effect is really due to the weak localization corrections and the present 
theory is not relevant to the experiment? 

However, it is not difficult to show that in the case under consideration the weak localization corrections originating 
from a contribution of Cooperons are totally suppressed by the magnetic field. This is not in contradiction with 
the fact that the system is close to the metal-insulator transition because strong localization is possible even if the 
Cooperons are absent. 

To calculate the contribution awL coming from the Cooperons we extended the standard derivation of the 
correctiorO to the case of the granulated metal. Using approximations developed previously one can obtain without 
difficulties the following expression for a 3-dimensional cubic lattice of the grains 



5a 



WL 



3 ^ 



■"/' Cq (0) , d3q 

cos qid- 



2'KVf) 



{2^r 



(10.3) 



where the function Cq (0) is the Cooperon taken at the frequency to — Q and quasi-momentum q. 

What remains to do is to take the e xpli cit expression for the Cooperon and compute the integral over t he q uasi- 
momentum q. However, if we use Eq. ( |3.6[ ) derived previously we get ide ntically zero. This is because Eq. (3.6) was 
deriv ed w ithout taking into account tunneling. We can modify Eq. (3^) writing as in Sec. IX the Dyson equation, 
Eq. (9.3), with E represented in Fig. 9. As a result, we obtain 
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Cq {iuon) = Ztti^o kn| H J5^{1- COS qid) + So [H) 



(10.4) 



In the limit under consideration, Eq. ( 3.14 ), the second term in the brackets in Eq. ( 10.4 ) is much smaller that 
the third one and one can expand the function Cq (iij^n), con sider ing the second te rm a s a perturbation. Restricting 
ourselves by the first order, substituting the result into Eq. ( 10. 3| ) and using Eq. ( 2.10|) we write the final result for 
the correction 5awL m the form 



5owL 

CTO 



8J 

TT 



8n{H) 



(10.5) 



Eq. ( 10. 5| ) shows that the weak localization correction in the strong magnetic fields considered here is always small. 
In contrast, the correction to the conductivity coming from the DOS 
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obtained in Sec. HI for T — can be considerably larger. The ratio of these two corrections takes ai H ^ H^ the 
form 



5awL ^ 24 fJ5\ ^^_i /Ao 
5(TDos TT VAo/ \J5 



(10.7) 



A 



which must be small in the limit J <C Ao/5. 

Now, let us estimate the corrections Saoes a-nd Sa\Y]^ using the parameters of the experimenlU. For the typical 
diameter 120 ± 20A of Al grains studied iru the mean level spacing 6 is approximately S ~ IK. Using the critical 
temperature Tc ~ 1.6K for Al we obtain for the BCS gap in a singl e grain the following result Ao ~ 1.8Tc « 3K. 
Substituting the extracted values of the parameters into Eq. ( 10. 6[) we can estimate the maximal increase of the 
resistivity. As a result, we obtain ('5/o//9o)„iax ~ ^-^^ which is somewhat smaller but not far from the value [Sp/ p)^-^„ ~ 1 
observed experimentally. 

Although our theory gives smaller values of ('5/9//Oo)max than the experimental ones, the discrepancy cannot be 
attributed to the weak localization effects. Using the experimental values of J, S and Ap we find from Eq. (10.7) that 
6awL is 10 times smaller than Saoos- The value of the correction 6awL/cro, Eq- (10.5), near He equals 2.8 x 10^^. 
Strictly speaking, all calculations have been done under the assumption of a large J ^ 1, while experimentally this 
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parameter is not large. This means that our theory does not take into account all effects that might be relevant for 
the experimentu. Possibly, at such small values of J as one has in the experiment charging effects become important 
reducing additionally the density of states. However, study of the effects of the Coulomb interaction is beyond the 
scope of the present pap er. 

We see from Eq. (3.25) that, if the orbital mechanism of the destruction of the superconductivity is more important 
than the Zccman one, then in the region of ultra high magnetic fields H ^ He the correction to the resistivity d ecay s 
as Sp/po ^ H^^ . In the opposite limiting case, when the Zeeman splitting is more important, we see from Eq. ( 7.12| ) 
that 5p/ pq ^ H^^. This means that the correction to the classical conductivity is still sensitive to the magnetic field 
far away from He- As concerns the weak localization correction 5awL given by Eq. (10.5), it decays at large magnetic 
fields as H~'^ and is always small. 

The dependence of the conductivity a on the magnetic field H is completely different at temperatures jT — T^ <C Tj, 
because in this region all types of the corrections considered in the previous sections can play a role depending on 
the value of the magnetic field. At low magnetic fields the Aslamazov-Larkin 5a al and anomalous Maki-Thompson 
^"iviT corrections give the main contribution because they are most divergent near the critical point. At the same 
time, as the magnetic field grows they decay faster than the correction Saoos originating from the reduction of the 
density of states. At a certain magnetic field H* all the corrections can become of the same order of magnitude. In 
this region of the fields, one can expect a negative magnetoresistance, although it cannot be large. Comparing all the 
types of the corrections with each other we estimate the characteristic magnetic field as H* ^ H^J^^ (Ao/(5). The 
theory presented here was developed under the assumption J^^(Ao/(5) ^ 1 and, hence, the characteristic field is very 
large, H* > H^. 

Experimentally, the peak in the resistivity at temperatures close to T^ can be estimated as 5p/ p^^ ~ 3.5 and is much 
smaller than the peak at low temperatures. This correlates, at least qualitatively, with our results. 



XI. CONCLUSION 

In this paper we presented a detailed theory of the new niechanism of the negative raagnetoresistance in granular 
superconductors in a strong magnetic field suggested recentlyQ to explain the experimenllj. We considered the limit of 
a large conductance, thus neglecting localization effects. It has been demonstrated that even if the superconducting 
gap in each granule is destroyed by the magnetic field the virtual Cooper pairs can persist up to extremely strong 
magnetic fields. However, the contribution of the Cooper pairs to transport is proportional at low temperatures to T^ 
and vanishes in the limit T ^ 0. In contrast, they reduce the one-particle density of states in the grains even at T = 0, 
thus diminishing the macroscopic conductivity. The conductivity can reach its classical value only in extremely strong 
magnetic fields when all the virtual Cooper pairs do not exist anymore. This leads to the negative magnetoresistance. 

We analyze both the orbital and Zeeman mechanisms of the destruction of superconductivity as well as the limits 
of low temperatures and temperatures close to the critical temperature T^.. The results demonstrate that, at low 
temperatures T <^ Tc, there must be a pronounced peak in the dependence of the resistivity on the magnetic field. 
This peak should be much smaller in the region of temperatures T <S^ Tc because, in this region the superconducting 
fiuctuations can contribute to transport, thus diminishing the role of the reduction of the density of states. We were 
able also to consider different regions of the magnetic fields. 
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FIG. 10. Resistivity of the granulated superconductors as a function of the magnetic field at fixed temperatures in different 
regimes. 

Qualitatively, the results are summarized in Fig. 9, where typical curves for the low temperatures T ^ Tc and 
temperatures close to T^ are represented. Both the functions reach asymptotically the value of the classical resistivity 
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Rq only at extremely strong magnetic fields. The resistivity R at low temperatures grows monotonously when 
decreasing the magnetic field. The function does not have any singularity at the magnetic field He destroying the 
superconducting gap A in a single grain. The real transition into the superconducting state occurs at a lower field 
Hc2- This field, in the region of parameters involved, is close to the field He- The resistivity R{H) remains finite as 
H — > Hc2 but its derivative diverges resulting in the infinite slope at H ~ Hc2 ■ The dependence of the resistivity R on 
the magnetic field at temperatures near Tc is more complicated. Already far from the field He the superconducting 
fluctuations start contributing to transport and the resistivity goes down when decreasing the magnetic field. A 
negative magnetoresistance is possible in this region only at magnetic fields H*, that can be much larger than the 
field He. 

The theory developed gives a good description of existing experiments. Although the experimental systems are 
close to the metal-insulator transition and localization effects as well as Coulomb interaction can play an essential 
role, our theory, where all these effects were neglected, gives reasonable values of physical quantities and allows to 
reproduce the main features of experimental curves. 

It was important for our calculations that the dimensionless conductance J 3> 1 of the sample was limited also 
from above, such that the inequality J ^ Ao/(5 was fulfilled. This means that the granulated structure of the 
superconductor was essential for us. However, the fact that the contribution of the superconducting fluctuations 
vanishes in the limit T -^ (AL and MT corrections are proportional to T^) seems to be rather general and not 
restricted by this inequality. Apparently, the negative magnetoresistance above the critical magnetic fleld persists and 
can be possible even in conventional bulk superconductors. We leave the region of large conductances J > Aq/S for 
a future study. 
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